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Abstract 

Early detection and accurate estimation of aortic stenosis (AS) severity are the most important predictors of successful long- 
term outcomes in patients. Current clinical parameters used for evaluation of the AS severity have several limitations 
including flow dependency. Estimation of AS severity is specifically challenging in patients with low-flow and low 
transvalvular pressure gradient conditions. A proper diagnosis in these patients needs a comprehensive evaluation of the 
left ventricle (LV) hemodynamic loads. This study has two objectives: (1 ) developing a lumped-parameter model to describe 
the ventricular-valvular-arterial interaction and to estimate the LV stroke work (SW); (2) introducing and validating a new 
index, the normalized stroke work (N-SW), to assess the global hemodynamic load imposed on the LV. N-SW represents the 
global hemodynamic load that the LV faces for each unit volume of blood ejected. The model uses a limited number of 
parameters which all can be measured non-invasively using current clinical imaging modalities. The model was first 
validated by comparing its calculated flow waveforms with the ones measured using Cardiovascular Magnetic Resonance 
(CMR) in 49 patients and 8 controls. A very good correlation and concordance were found throughout the cycle (median 
root mean square: 12.21 mL/s) and between the peak values (r = 0.98; SEE = 0.001, p<0.001). The model was then used to 
determine SW using the parameters measured with transthoracic Doppler-echocardiography (TTE) and CMR. N-SW showed 
very good correlations with a previously-validated index of global hemodynamic load, the valvular arterial impedance (Z VA ), 
using data from both imaging modalities (TTE: r = 0.82, SEE = 0.01, p<0.001; CMR: r = 0.74, SEE = 0.01, p<0.001). 
Furthermore, unlike , N-SW was almost independent from variations in the flow rate. This study suggests that considering 
N-SW may provide incremental diagnostic and prognostic information, beyond what standard indices of stenosis severity 
and provide, particularly in patients with low LV outflow. 
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Introduction 

Aortic stenosis (AS) is a complex "systemic" disease. There are 
compelling epidemiological and histopathological data suggesting 
that "degenerative" calcified AS is, in fact, an active and 
multifaceted disease that involves atherosclerotic-like and elasto- 
calcinosis-like processes [1,2,3]. 

The degenerative process is not limited to the aortic valve; it 
involves the vascular system distal to the valve as well, which in 
turn can further contribute to deteriorating of left ventricle (LV) 
function [4]. 

Moreover, in patients with severe AS, "paradoxical" low flow 
(stroke volume index (SVi) <35 mL/m 2 ) and consequently low 
transvalvular pressure gradients (<40 mmHg) despite the presence 
of preserved LV ejection fractions (>50%), is a challenging clinical 
entity [5,6,7]. This mode of presentation of severe AS is relatively 



frequent (up to 35% of cases) and has been shown to be associated 
with a more advanced stage of the disease [5,6,7]. Yet, a majority 
of these patients do not undergo surgery due to inadequate 
diagnosis. In these patients, the presence of low pressure gradient 
in conjunction with a normal LV ejection fraction may easily lead 
to an underestimation of AS severity [5,6,7]. Prognosis in such 
patients is usually poor (survival rates <50% at 3-year follow-up) if 
treated medically and operative risk is high (up to 33%) if treated 
surgically [5,8]. Early detection and accurate estimation of AS 
severity are therefore of primary importance. For proper diagnosis 
of these patients, a comprehensive evaluation of LV hemodynamic 
loads is crucial. Such a comprehensive evaluation helps to identify 
needs for closer follow-ups, and to improve risk stratification and 
clinical decision making [8]. 
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Figure 1. Schematic diagrams, (a) electrical representation, (b) schematic representation of the lumped parameter model used to simulate left- 
sided heart in presence of aortic stenosis and/or systemic arterial hypertension. LV: left ventricle, AS: aortic stenosis, E| v (t): normalized time-varying 
elastance (please see Table 1 for all other parameters used in the lumped parameter model). 
doi:10.1371/journal.pone.0086793.g001 



Recently, Briand et al. (2005) proposed the valvulo-arterial 
impedance (Zva) (equation 8) which represents the valvular and 
arterial factors that oppose ventricular ejection by absorbing the 
mechanical energy developed by the LV. This index has been 
shown to be a strong predictor of clinical outcomes in 
asymptomatic AS patients and is superior to the standard indices 
of AS severity in predicting LV dysfunctions [9, 1 0] . However, Zva 
is flow-dependent and therefore has limited use specifically in AS 
patients with low flow and low pressure gradient conditions [8] . 

The LV stroke work has been shown to be effective in 
characterizing the LV loads and consequently in characterizing 
patient's outcome and in assessing the inotropic state in patients 
with AS [11,12,13]. LV stroke work is the energy that the ventricle 
delivers to the blood at ejection, and potential energy, necessary to 
overcome the viscoelastic properties of the myocardium itself. 
Therefore, A method that would allow simple, non-invasive, and 
accurate estimation of the LV stroke work would contribute 
towards improving the clinical management of patients with AS 



and especially the ones with paradoxical low flow. However, 
estimation of this parameter requires the knowledge of instanta- 
neous LV pressure and volume. Cardiac catheterization with 
conductance catheters can provide this information but this 
method is invasive, expensive and may cause cerebral embolisms 
[14], A non-invasive alternative to estimate LV stroke work is to 
model the cardiovascular system using a lumped parameter model 
based on non-invasive recorded inputs. 

The first objective of this study was therefore to develop and to 
introduce a simple lumped parameter model, solely based on non- 
invasive parameters, to describe the ventricular-valvular-arterial 
coupling and to investigate the impact of AS and vascular diseases 
on LV workload. The second objective was to introduce and 
validate a new index easily measurable in the clinical setting, i.e., 
the normalized LV stroke work (N-SW). This index in J/mL, 
representing the energy required to eject 1 mL of blood through 
the valvulo-arterial system, is calculated by dividing the estimated 
LV stroke work by the stroke volume (SV). 
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Table 1. Summarized cardiovascular parameters used to simulate all cases. 





Description 


Abbreviation 


Value 


Aortic valve parameters 


Effective orifice area 


EOA 


From TTE and CMR data 


Aortic cross sectional area 




From TTE and CMR data 


Valvular energy loss coefficient 


E L Co 


(EOA)A 
A -EOA 


Variable aortic valve resistance 


Rav 


IElCo 1 6 


Aortic valve inductance 


Lav 


VEZCo 


Systematic circulation parameters 


Aortic resistance 


Rao 


0.05 mmHg.s.mr 1 


Aortic compliance 


Cao 


0.5 ml/mmHg 


Systemic vein resistance 


Rsv 


0.05 mmHg.s.mr 1 


Systemic arteries and veins compliance 


CsAC 


Initial value: 2 ml/mmHg 

Adjust for each degree of hypertension 


systemic arteries resistance 

(including arteries, arterioles and capillaries) 




Initial value: 0.8 mmHg.s.mr 1 

Adjust according to the obtained total systemic resistance 


Output condition 


Central venous pressure 


Pcvo 


4 mmHg 


Input condition 


Mitral valve mean flow rate 


Qm„ 




Other 


Constant blood density 




1 050 kg/m 3 


Cardiac output 


CO 


From TTE and CMR data 


Heart rate 


HR 


From TTE and CMR data 


Duration of cardiac cycle 


T 


From TTE and CMR data 



doi:1 0.1 371 /journal.pone.0086793.t001 



Methods 

Lumped Parameter Modeling 

A schematic diagram of the lumped parameter model is 
presented in Fig. 1 [15]. This model includes three different sub- 
models: 1) LV model; 2) AS model; 3) systemic circulation model. 
The effect of AS on LV workload was investigated under several 
numerical conditions using datasets obtained non-invasively by 
transthoracic echocardiography (TTE) and cardiovascular mag- 
netic resonance (CMR) in healthy subjects and in AS patients. All 
other parameters used in the lumped parameter model are listed in 
Table 1. 

Heart-arterial model. Heart function was described by time 
varying elastance as following 



E(ty- 



P L v(t) 
V(t)-V 0 



(1) 



where PlvW, V(t) and V 0 are the left ventricular pressure, the 
left ventricular volume and the unloaded volume [15], 
respectively. The amplitude of E(t) was normalized with respect 
to maximal elastance E max , i.e., the slope of the end-systolic 
pressure-volume relation, giving E N (t N ) = E(t)/E max . Time then 
was normalized with respect to the time to attain peak 
elastance, T Emax (t N = t/T Emax ). Normalized time-varying ela- 
stance curves En^n) have similar shapes in the normal human 
hearts with various inotropic situations or for diseased human 
hearts despite the existence of differences with regard to etiology 



of cardiovascular diseases 
elsewhere [15] 



[16,17]. More details can be found 



E m axEN(t/T E 



V(t)-V 0 



(2) 



Modeling aortic stenosis. Aortic stenosis was modeled using 
the semi-analytical formulation for the net pressure gradient 
(TPG„ et ) across the stenotic valve during LV ejection introduced 
by Garcia et al. [18]. This formulation expresses the instantaneous 
net pressure gradient across the stenotic valve (after pressure 
recovery) as a function of the instantaneous flow rate and the 
energy loss coefficient and links the LV pressure to the ascending 
aorta pressure (P 'lv(I) — P U(0) 



TPG na =P L v{i)-P A {t)-- 



and 



2np dQ(t) p 
s/E^Co dt 2E L Co 2 



E L Co-- 



(EOA)A 
A -EOA 



Q 2 ^) (3) 



(4) 



where E^Co, EOA, A, p and Q are the valvular energy loss 
coefficient, the effective orifice area, the aortic cross sectional area, 
the fluid density and the transvalvular flow rate, respectively. 
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Table 2. Baseline characteristic data. 







Healthy Subjects 


AS Patients 




(n = 8, mean ± SD) 


(n = 49, mean ± SD) 


Patient description 


Age (years) 


34±8 


63±16* 


Sex (Men %) 


75 


63 


Body surface area (m 2 ) 


1.93±0.26 


1.82±0.19 


Left ventricle function and geometry 


Left ventricle end-diastolic internal dimension (mm) 


50±3 


45±5 


Interventricular septal thickness (mm) 


9±1 


12±2* 


Left ventricle mass (g) 


230±54 


194±55 


Left ventricle mass index (g/m 27 ) 


52±14 


49±52 


Left ventricle ejection fraction (%) 


66±4 


66±5 


Arterial hemodynamics 


Systemic arterial compliance (mL.m _2 .mmHg _1 ) 


1.17±0.32 


0.78 ±0.26* 


Systemic vascular resistance (dyne.s.cm -5 ) 


1335 ±284 


1801 ±774 


Systolic arterial pressure (mmHg) 


116±10 


128±22 


Diastolic arterial pressure (mmHg) 


77±5 


71 ±11 


Mean arterial pressure (mmHg) 


90±7 


90±13 


Valvulo-arterial impedance (mmHg/mL/m 2 ) 


3.35±0.71 


3.66±0.85 


Valve hemodynamics 


Mean transvalvular gradient (mmHg) 


5±1 


21±11* 


Effective orifice area (cm 2 ) 


2.67±0.47 


1.31 ±0.59* 


Aortic diameter (mm) 


31 ±4 


32±4 


Left ventricular stroke work 


Stroke work (J) 


0.85 ±0.2 


1.45±0.75 


Normalized stroke work (J/mL) 


0.009±0.002 


0.01 7± 0.007 



*p<0.05 with healthy. 

doi:1 0.1 371 /journal.pone.0086793.t002 



Computational algorithm. The lumped model illustrated in 
figure 1 was analyzed numerically by creating and solving a system 
of ordinary differential equations in Matlab Simscape (Math- 
Works, Inc). Capabilities of this program were enhanced by 
adding additional codes to meet demands of cardiac circuit. 
Fourier series representation of experimental normalized elastance 
curve for human adults was used to generate a signal to be fed into 
the main program. Equation 3, representing the transvalvular 
pressure gradient across the aortic valve, was represented by an 
inductance and a variable resistor as depicted in figure 1. 
Simulation started at the onset of isovolumic contraction and the 
elastance signal drives the program by feeding elastance value 
related to each time step in the cycle to the equation 1. The left 
ventricle volume V(t) was calculated using the left ventricle 
pressure Plv an d elastance values by equation 1. The Plv used at 
the beginning of calculation was the initial value assumed across 
the variable capacitor and was automatically adjusted later by the 
system of equations as solution advances. The left ventricle flow 
rate subsequendy was calculated as the time derivative of the left 
ventricle volume. After few initial cycles, solution converged. A 
diode with very low on resistance and off conductance was used in 
the aortic valve to prevent backflow from the valve. Matlab's 
"ode23t" trapezoidal rule variable-step solver was used to solve 
system of differential equations with initial time step of 0.1 
milliseconds. The Convergence residual criterion was set to 10 



and initial voltages and currents of capacitors and inductors set to 
zero. 

Determining arterial compliance and peripheral 
resistance. The total systemic resistance was computed as the 
average brachial pressure over the cardiac output (assuming a 
negligible peripheral venous pressure (mean ~5 mmHg) com- 
pared to aortic pressure (mean — 100 mmHg)). This total systemic 
resistance represents the electrical equivalent resistance for all 
resistances in the current model. Because what the left ventricle 
faces is the total systemic resistance and not the individual 
resistances, for the sake of simplicity, we elected then to consider 
the aortic resistance, R ao , and systemic vein resistance, Rsv, as 
constants and to adjust the systemic artery resistance, Rsa, 
according to the obtained total systemic resistance. 

Physiologically, arterial hypertension is determined by two 
factors [19]: (1) a reduction in the caliber of small arteries or 
arterioles with an ensuing increase in systemic vascular resistance 
and mean blood pressure, and (2) a reduction in the arterial 
compliance with a resulting increase in pulse pressure (systolic 
minus diastolic blood pressure). In this study, we fitted the 
predicted pulse pressure to the actual pulse pressure (known by 
arm cuff sphygmomanometer) by adjusting the systemic compli- 
ance (Csac)- Therefore, C SAC adjustment was done by a simple 
trial and error for each degree of hypertension [20] . 
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Figure 2. Schematic diagram. LV stroke work and N-SW. 
doi:1 0.1 371 /journal.pone.0086793.g002 



In vivo Measurements 

Ethics statement. This study was approved by the Ethics 
Committee of Research of Laval Hospital affiliated to Laval 
University. Informed consents were obtained from all patients. 

Study population. Forty-nine patients with mild to severe AS 
(63% men, age 63± 16 years) and eight healthy subjects (75% 
men, age 34 ±8 years were included in this study (Table 2). 
Exclusion criteria were: (1) age <21 years old. Patients with AS are 
often older than 21 years, also younger patients are treated as 
pediatric population according to the American Heart Association 
guideline; (2) LV ejection fraction <50%. SV measurements in 
patients with reduced LV ejection fraction may show inconsisten- 
cies; (3) moderate or severe mitral or aortic regurgitation. SV 
measurements in these patients may show inconsistencies; (4) poor 
TTE data quality and standard contra-indications to CMR data. 
To demonstrate the feasibility of our method, only patients with 
useable data were considered in our study. In 33% of AS patients, 
the valve morphology was bicuspid. Other measured character- 
istics of patients were: height (166±9 cm), weight (76±13 kg), 
waist circumference (95 ±11 cm), blood pressure 
(128±22 mmHg/71±ll mmHg), body surface area 
(1.82±0.23 m 2 ), diabetes (16%) and metabolic syndrome (35%). 
The summary of performed measurements is presented in Table 2. 
All patients provided written informed consent under the 
supervision of the Institutional Review Board. The initial AS 
severity classification at study entry was based on TTE-derived 
effective orifice area (EOA): normal (EOA>2.0 cm 2 ), mild 



(1.5<EOA<2 cm 2 ), moderate (1.0 cm 2 <EOA<1.5 cm 2 ) and 
severe (EOA<1.0 cm 2 ). 

Transthoracic Doppler-echocardiography (TTE). TTE 

exams were performed and analyzed by two experienced 
echocardiographers, and conducted according to the American 
Society of Echocardiography guidelines [21]. They included: 

(a) Valve hemodynamic parameters: the left ventricle outflow 
track (LVOT) diameter, LVOT flow velocity measured by 
pulsed-wave Doppler, the aortic transvalvular jet velocity 
measured by continuous-wave Doppler and valve effective 
orifice area (EOA) using the continuity equation as follow, 



(b) 



EOAj 



sv Ll 



VTI A o 



( VTIlvot x Alvot) 
VTI AO 



(5) 



where SVlvot, Alvot and VTIlvot are the stroke 
volume measured in the LVOT, the cross-sectional area 
of the LVOT calculated assuming a circular shape 
(LVOT diaml , tel . 2 x0.785), and the velocity-time integral of 
the LVOT, respectively. 

Parameters of LV systolic function: the left ventricle ejection 
fraction (LVEF), Sa wave and left ventricle geometry [22]; 
Parameters of LV diastolic function: Ea wave, El A ratio, left 
ventricle geometry and mass indexed to height 2 ' 7 (LVMi) 
[23,24]; 
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Figure 3. Comparison of simulated and CMR flow waveforms, (a) healthy subject, (b) patient with mild AS, (c) patient with moderate AS, (d) 

patient with sever AS. 

doi:1 0.1 371 /journal.pone.0086793.g003 



(d) Vascular hemodynamic parameters: the systemic arterial 
compliance (SAC) and the systemic vascular resistance (SVR) 

[9];. 



LVSP _ P S + (TPG) A 

~sv~~ 



SVi 



(8) 



SAC=SVi/PP 
SVR = WMAP/CO 



(6) 
(7) 



where SVj, PP, A-LAP, and CO are stroke volume indexed by the 
body surface area, pulse pressure, mean arterial pressure and 
cardiac output, respectively. 

(e) Global (valvulo+arterial) LV load, estimated by the valvulo- 
arterial impedance (Z V a) [9], as follow: 



where LVSP, P$ and (TPG) AS are LV systolic pressure, 
aorta systolic pressure and mean transvatvular pressure 
gradient, respectively. 

Cardiovascular Magnetic Resonance (CMR). The CMR 

study was performed 2 to 4 weeks after the TTE study with the use 
of a 1.5 Tesla scanner and a dedicated cardiac phase-array 
receiver coil (Achieva, Philips Medical Systems, Best, The 
Netherlands). CMR image acquisitions and analyses were 
performed by investigators blinded to clinical and TTE results, 
as previously described [25,26,27]. 

A custom-made research application was developed using 
Madab software (Mathworks, Natick, Ma) to process and analyze 
CMR images [24,25,26]. The valve EOA from CMR (EOAqmr) 
was then calculated as follow: 
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Figure 4. Correlation between results based on TTE and CMR measurements, (a) correlation between computed normalized left ventricular 
stroke works based on parameters measured by TTE and CMR, (b) correlation between computed valvulo-arterial impedance based on parameters 
measured by TTE and CMR. 
doi:1 0.1 371 /journal.pone.0086793.g004 



EOA, 



SV C 



VTI A 



(9) 



where SVcmr is the stroke volume using 1/3 Simpson's rule to 
integrate the systolic flow and VTIao is the velocity-time integral 
of the peak aortic flow velocity measured at 1 0 mm downstream of 
the valve during systole [25,27,28,29]. 

Computed normalized stroke work index (N-SW). We 
have introduced a lumped-parameter method for clinical practice 
that accurately investigates the impact of AS and concomitant 
vascular diseases on the LV function. This method only needs few 
non-invasively measured quantities described as follows. TTE and 
CMR can both provide the effective orifice area [El Co) and aortic 
cross sectional area (A). Equation 4 uses these quantities to 
calculate the valvular energy loss coefficient (ElCo). Subsequendy 



Rav £ 



2E,Co 2 ' 



and L av ( 



2np 
\/E L Cd 



in Equation 3 and Fig. 1 are 



calculated. The total systemic resistance is computed as the 
average brachial pressure over the cardiac output and is used to 
calculate RgA assuming constant values for other resistances in the 
circuit. The systemic compliance (Csac) was adjusted by fitting the 
predicted pulse pressure to the measured pulse pressure. The LV 
stroke work, representing the work of the left ventricle during each 
heart beat (J Pz,vdVi,y), was then computed using the left 
ventricle volume, V(t), and the left ventricle pressure, Plv- 
Normalized LV stroke work (N-SW) is calculated by dividing the 
estimated LV stroke work by the SV (see Figure 2 for schematic 
diagram). 

Statistical analysis. Results were expressed as mean ± SD. 
CMR and TTE measurements were compared by 2-tailed paired 
Student t-tests or one-way ANOVA when appropriate. Association 
and agreement between variables were assessed by Pearson's 
correlation and Bland-Altman analyses, respectively. Statistical 
analysis was performed with SPSS 1 7 (SPSS, Chicago, IL). 
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Figure 5. Correlation between normalized left ventricular stroke work and valvulo-arterial impedance, both computed using the 
lumped parameter model based on the parameters measured by (a) TTE and (b) CMR. 
doi:1 0.1 371 /journal.pone.0086793.g005 
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Figure 6. Dependency of valvulo-arterial impedance to effective orifice area and stroke volume, (a) as a function of effective orifice area 

(EOA), (b) as a function of stroke volume for normal arterial pressure. 

doi:10.1371/journal.pone.0086793.g006 



Results 

Validation of the Lumped Parameter Model using in vivo 
CMR Flow Wave Forms 

Using the above mentioned protocol, data obtained non- 
invasively in forty-nine patients with mild to severe AS and eight 
healthy subjects (Table 2) were incorporated in the model and the 
LVOT flow waveforms resulted from the lumped parameter 
model were compared with in vivo CMR flow waveforms with an 
average root mean square error of 12.21 mL/s. There was a very 
good correlation between the peak values of flow waveforms 
obtained from simulations and CMR measurements (r = 0.98; 
SEE = 0.001, p<0.001). Figure 3 presents few examples of 
comparison between the simulated and CMR flow waveforms. 



Estimation of Normalized Stroke Work Based on the 
Input Data from TTE and CMR 

Figure 4 shows a strong correlation between the results 
calculated using lumped parameter model based on TTE and 
CMR imaging datasets. Fig. 4a shows the correlation for the 
computed N-SW (r = 0.9, SEE = 0.005, p<0.001) and Fig. 4b 
shows the correlation for the Zva (r = 0.78, SEE = 0.05, p<0.001). 

Correlation Analysis of the Normalized Stroke Work 

Figure 5 shows the strong correlation (TTE, Fig. 5a: r = 0.82, 
SEE = 0.01, p<0.001; CMR, Fig. 5b: r = 0.74, SEE = 0.01, 
p<0.001) between N-SW and 7.va, which has been shown to be 
a powerful predictor of mortality in patients with AS. 




SV: 30 mL 
SV: 50 mL 
SV: 70 mL 
SV: 90 mL 



0.5 



£ 



o 



1.5 2 2.! 

Effective orifice area (cm 2 ) 



0.03 



0.025 



0.02 



0.015 



0.01 



0.005 



20 




-«-EOA= 2 cm2 
^-EOA=1.5 cm2 
-e-EOA=1 cm2 
-&-EOA=0.5 cm2 



-»-EOA=1.7 cm2 
-»-EOA=1.2 cm2 
-»-EOA=0.7 cm2 
-±-EOA=0.3 cm2 



40 



60 



80 100 
Stroke Volume (mL) 



Very severe Severe Moderate Mild 



Figure 7. Dependency of normalized left ventricular stroke work to effective orifice area and stroke volume, (a) as a function of 
effective orifice area (EOA), (b) as a function of stroke volume for normal arterial pressure. 
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Figure 8. Valvulo-arterial impedance and normalized left ventricular stroke work, (a) Valvulo-arterial impedance (Z VA ) as a function of 
stroke volume, (b) normalized left ventricular stroke work (N-SW) as a function of stroke volume. In both panels a subset of patients with low flow and 
low pressure gradient conditions (Table 2) were considered. In these patients, AS was either moderate (EOA= 1.48 cm 2 , n = 4) or severe 
(EOA = 0.85 cm , n = 5). In both groups, the stroke volume index (SVi) was less than 35 ml/m 2 and the transvalvular pressure gradient was less than 
20 mmHg. Both Z VA and N-SW were computed using the lumped parameter model based on the parameters measured by CMR. 
doi:1 0.1 371 /journal.pone.0086793.g008 



Effect of the Flow Rate Condition on Normalized Stroke 
Work and Valvulo-arterial Impedance Based on TTE and 
CMR Input Data 

Figures 6 and 7 show the effect of the flow rate condition on 
ZyA and N-SW. These figures illustrate hypothetical cases for a 
large range of stroke volumes: from 30 and 90 mL and with 
different AS severities: EOA of 0.3 to 2 cm 2 . Figure 6 demon- 
strates that Zyyj is significantly influenced by variations in flow 
rates: ZyA decreased when the stroke volume was reduced from 90 
to 30 mL. In contrast, as Figure 7 shows, the flow rate had 
minimal effects on N-SW. 

Discussion 

Advantages of Normalized Stroke Work 

Patients with AS often have concomitant arterial diseases. This 
should alert clinicians that the degenerative process is not only 
limited to the aortic valve, but also involves the arterial system 
distal to the valve. The latter can further contribute to 
deteriorating LV function [2,30]. Under this condition, the LV 
faces an increased double load: a valvular load imposed by the AS 
plus an arterial load. Therefore, it is important to assess the 
hemodynamic load imposed on the LV. To this effect, ZyA 
representing the valvular and arterial impedances opposing 
ventricular ejection, was recently proposed [4]. 

N-SW represents the global hemodynamic load that the LV 
faces for ejecting a unit volume of the blood. It has a very good 
correlation with ZyA . However, using N-SW to assess the AS has 
privilege to ZyA for two main reasons: 1) N-SW is not flow 
dependent while, by nature, Zva is. This may limit the utility of 
ZyA in patients with low-flow and low-gradient AS [8], 2) N-SW 
determines the actual total mechanical load imposed to the LV, 
while ZyA gives an estimate of that load. 



Estimation of Normalized Stroke Work and Valvulo- 
arterial Impedance by TTE versus CMR 

This study demonstrates a very good correlation between the N- 
SWs computed based on the data from TTE and CMR methods. 
Similarly Zy A computed based on the data from TTE and CMR 
methods exhibit a very good correlation. TTE is the primary 
imaging technique for the assessment of AS severity and its 
progression. When TTE measurements are inconclusive or show 
discordances, AS severity should be confirmed with other 
techniques either invasively using cardiac catheterization or non- 
invasively using CMR [26,31]. Results of this study suggest that 
there is a good agreement between TTE and CMR for the 
calculation of N-SW. 

Clinical Implications 

Since the management of patients with AS remains a source of 
debate [3,5,6,7], supplementing currently available indices with 
the N-SW could be useful in identifying patients who are at high 
risks and are in need of closer follow-ups and further investiga- 
tions. 

More specifically, as already reported in [6], a significant 
number of patients with severe AS have low transvalvular flow 
rates (SVi <35 ml/m 2 ), low transvalvular pressure gradients 
(<40 mmHg) and preserved LV ejection fractions (>50%). 
Clinically, this highly insidious situation represents a challenging 
clinical entity because the AS may appear less severe on the basis 
of low transvalvular pressure gradients. This situation may lead 
clinicians to erroneously conclude that no surgery is required 
because the stenosis is not considered to be severe. Indeed, these 
patients are at a more advanced stage of the disease. When 
compared with patients with severe AS but normal LV output, 
thus high pressure gradients, these patients are characterized by 
higher degrees of LV concentric remodeling, lower LVEFs and 
reduced mid-wall shortenings. It was found that these patients 
have poor prognoses and outcomes [6,7]. Therefore, it can be 
conceived that a longstanding increase in the LV load may result 
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in a more pronounced LV concentric remodeling, a smaller LV 
cavity size, and a decrease in intrinsic myocardial function. Early 
detection and accurate estimation of AS severity are, therefore, of 
primary importance. In this regard, for a comprehensive 
evaluation of the valvular and vascular loads, calculation of 2va 
appears to be particularly useful for patients with AS. However, in 
low-flow low-gradient AS patients, Zva would likely be less precise 
in the assessment of the LV global hemodynamic load, due to its 
high flow-dependency (Figure 8a) [8]. In these patients, N-SW 
which is less flow-dependent, may be useful to better quantify the 
LV hemodynamic load and to better identify the actual disease 
severity that is masked by this paradoxical low flow phenomenon 
(Figure 8b). 

Conclusions 

Patients with AS often have concomitant arterial diseases. 
Under this condition, the LV faces a double load: a valvular 
load imposed by the AS plus an arterial load. The double load 
results in increasing of the LV stroke work. We, therefore, 
introduced a lumped parameter method, capable of accurately 
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